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O Abstract 

An application of the thermodynamic field theory (TFT) to trans- 
<— | port processes in L-mode tokamak plasmas is presented. The nonlin- 

Q_i ear corrections to the linear ("Onsager") transport coefficients in the 

collisional regimes are derived. A quite encouraging result is the ap- 
pearance of an asymmetry between the Pfirsch-Schluter (P-S) ion and 
^ electron transport coefficients: the latter presents a nonlinear correc- 

tion, which is absent for the ions, and makes the radial electron coeffi- 
cients much larger than the former. Explicit calculations and compar- 
isons between the neoclassical results and the TFT predictions for JET 
plasmas are also reported. We found that the nonlinear electron P-S 
transport coefficients exceed the values provided by neoclassical theory 
by a factor, which may be of the order 10 2 . The nonlinear classical 
coefficients exceed the neoclassical ones by a factor, which may be of 
order 2. The expressions of the ion transport coefficients, determined 
. — , by the neoclassical theory in these two regimes, remain unaltered. 

The low-collisional regimes i.e., the plateau and the banana regimes, 
are analyzed in the second part of this work. 

1 Introduction 

The thermodynamic field theory (TFT) was proposed in 1999, to describe 
the behaviour of thermodynamic systems beyond the linear (" Onsager" ) re- 
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gion [1!. Attempts to derive a generally covariant thermodynamic field the- 
ory (GTFT) can be found in refs |2j. The TFT extends the theory previously 
formulated by Prigogine in 1954, which was applied only to thermodynamic 
systems close to equilibrium. The characteristic feature of this theory is 
its purely macroscopic nature. We do not mean a formulation based on the 
macroscopic evolution equations, but rather a purely thermodynamic formu- 
lation starting solely from the entropy production and from the transport 
equations, i.e., the flux-force relations. The latter provide the possibility of 
defining an abstract space, (the thermodynamic space) whose metric is given 
by the transport matrix. The law of evolution is not the dynamical law of 
motion of particle, or the set of two-fluid macroscopic equations of plasma 
dynamics. The evolution in the thermodynamic configurations is rather de- 
termined by postulating three purely geometrical principles: the Shortest 
Path Principle, the Closeness of the Thermodynamic Field Strength, and 
the Principle of Least Action. From theses principles, a set of field equa- 
tions, constraints, and boundary conditions are derived. These equations, 
referred to as the thermodynamic field equations, determine the nonlinear 
corrections to the linear ("Onsager") transport coefficients. 
The validity of this theory in the weak-field approximation has been suc- 
cessfully tested in many examples of a — a and (3 — (3 processes, such as the 
thermoelectric effect and the unimolecular triangular chemical reaction [Tj^j 
The thermodynamic field equations, in the weak-field approximation, have 
also been applied to several a — (3 processes. For example, the Field-Koros- 
Noyes model, in which the thermodynamic forces and flows are related by an 
asymmetric tensor, was analyzed in ref. |3]. Even on this case, the numerical 
solutions of the model are in agreement with the theoretical predictions of 
the TFT. More recently, the Hall effect [5] has been analyzed in the nonlin- 
ear region. In each of these papers, it was shown that the TFT successfully 
describes the known physics in the nonlinear region and also predicts new 
interesting effects, such as the nonlinear Hall effect. The theoretical predic- 
tions of the nonlinear Hall effect have been confirmed experimentally [6] . 
The TFT has also been used to study transport processes in magnetically 
confined plasmas. Preliminary results can be found in refs |7j and [5]. The 
study of the behaviour of a plasma in the presence of an inhomogeneous and 
curved magnetic field is one of the main objects of the neoclassical theory. 

1 Here, we adopt the terminology of De Groot and Mazur [3], i.e., when the velocity 
distribution function is an even (odd) function of the velocities of the particles, a processes 
is said to be an a-processes (/3-processes). It is possible to show that this definition 
implies that a processes only involve the symmetric part of the Onsager tensor, whereas 
f3 processes only involve the skew-symmetric one. 
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One of the most important results of the neoclassical theory is that the global 
geometry of the magnetic field has a very strong influence on the transport 
processes. The neoclassical theory is able to derive the transport coefficients 
for ions and electrons when the plasma is magnetically confined in tokamak 
reactors. In spite of its elegant and coherent formulation, the theoretical 
predictions of the neoclassical theory are in strong disagreement with ex- 
perience. The experimental ion heat flux measured in tokamak plasmas is 
roughly in agreement with the neoclassical theory. However, the electron 
particle flux and the electron heat flux are about 10 2 -j- 10 3 greater than 
the values computed by the neoclassical theory. This difference between the 
experimental and the neoclassical flux is referred to as the anomalous flux. 
For many physicists, the origin of this discrepancy is mainly attributed to 
turbulent phenomena existing in tokamak plasmas. Fluctuations in plasmas 
can become unstable and therefore amplified. According to this interpreta- 
tion, fluctuations will successively interact in a nonlinear way leading the 
plasma to a state, which is far away from equilibrium. In this condition, the 
transport properties are supposed to change significantly. 
In this one and the successive paper, we shall analyze the influence of the 
nonlinear contributions, evaluated by the thermodynamic field theory, on 
the transport processes in tokamak plasmas. As mentioned above, the neo- 
classical theory (i.e., the linear theory) fails with a factor 10 2 -j- 10 3 and 
thus magnetically confined plasmas represent an ideal case for testing the 
validity of our thermodynamic approach. For simplicity, in our calculations 
we deal with fully ionized plasmas defined as a collection of electrons and 
positively charged ions |10| . In the present and in the subsequent work j9], 
we analyzed in great detail transport processes in tokamak plasmas in the 
collisional and in low-collisional regimes. In this first part, we limit our- 
selves to a confined plasma in collisional regimes, i.e., in the classical and 
Pfirsch-Schluter regimes, providing the complete set of the ion and electron 
nonlinear transport equations and analyzing the solution for JET-plasmas. 
We report a new set of nonlinear transport equations for plasmas in colli- 
sional regimes different from those established in refs [7] and |8]. This is 
due to the fact that here we adopt a definition of the Pfirsch-Schluter ther- 
modynamic forces different from that reported in refs |7j and [8] . The work 
is organized as follows: in sections ([2]) and Q we derive the expressions 
of the entropy productions for ions and electrons in magnetically confined 
plasmas and the transport equations in the linear region. After a brief sum- 
mary of the general theory of the thermodynamic field (TFT), we adapt the 
formalism in order to obtain the thermodynamic field equations for a fully 
ionized plasma. The solution of these equations determine the nonlinear 
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corrections to the linear ("Onsager") transport coefficients. The analysis of 
the classical and Pffrsch-Schliiter transport is conducted in sections (|5| and 
in section ^ we derive the gauge-invariant form of the solutions, which 
are explicitly analyzed in section Q. As general result we find that, in the 
Pfirsch-Schluter regime, the thermodynamic fluxes are linked to the thermo- 
dynamic forces by an amplification factor J- ps times the Onsager matrix. In 
order to interpret the mathematical results found in section Q in terms of 
collisional mechanisms, a kinetic model can be used originally introduced in 
ref. [TT]. Through this model, it was possible to show that, in the collisional 
regimes, the transport coefficients are approximatively given by the linear 
(Onsager) transport coefficients times a function, which is proportional to 
the inverse of the electron collision time r e . At the end of section ^ we 
also find specific calculations for JET plasmas. The main conclusion of our 
analysis is that the electron nonlinear Pfirsch-Schluter transport coefficients 
exceed the values provided by neoclassical theory by a factor J- ps , which may 
be of order 10 2 . The values of the ion transport coefficients remain, how- 
ever, unaltered. These results are in line with experimental observations. 
The definition of gauge invariance of the field equations for our problem 
and the proof that our choice of boundary conditions respects the princi- 
ple of covariance is reported in appendix Details on the derivation of 
the appropriate boundary conditions for our problem can be found in ap- 
pendix (|9]). The explicit solution of the field equations, submitted to the 
appropriate boundary conditions, is given in appendix (9.1) 



2 Entropy Production of Magnetically Confined 
Plasmas 

We consider a plasma consisting of electrons and a single species of ions, in 
the presence of an external axisymmetric confining magnetic field. As it is 
known, neglecting the drift and modified drift mechanism, in magnetically 
confined plasmas, we have three distinct mechanisms of transport: banana, 
plateau and the Pfirsch-Schluter transport regimes (see, for example, ref. 
|12j). In ref. |12] we can find the evaluation of the entropy production a a 
of species a (a = e for electrons and a = i for ions). Let us first introduce 
the dimensionless (density of) entropy productions of species a, S a : 

£<* = ^ a a (1) 
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where T a and n a indicate the relaxation time and the number density of 
particles of species a, respectively. The expression for T, a is [12] : 

N 

E Q = -g«T Q G?« - ^ <£( 2n+1 >T a G?( 2n+1 > (2) 
n=l 

The upper limit is oo, but in practice the sum is truncated at n = N. 
The Gr ( 2n+1 ) and the g"^ 2n+1 ^ indicate the dimensionless collisional terms 
(often also referred to as the dimensionless generalized frictions) and the 
dimensionless irreducible Hermitian moments, respectively, is not an 
Hermitian moment (i.e., a moment defined by the deviation from the local 
Maxwellian) but indicates the dimensionless electric current density. Index 
r denotes the components of this vector. Expression ^ is completely gen- 
eral: it is based only on the assumption that the deviation from the local 
equilibrium is small and that tensor Hermitian moments provide negligi- 
ble contributions. The specificity of the neoclassical theory appears when 
we express the generalized frictions Gr ( 2n+1 * ) in terms of the dimensionless 
Hermitian moments qy^ 2n+1 ^ by making use of the approximated moment 
equations. These equations are obtained assuming the validity of the drift 
approximation. Indicating with e the drift parameter i.e., the ratio of the 
Larmor radius to the length of the gradient of the magnetic field, all terms 
of order e 2 or higher appearing in the moment equations can be neglected. 
This means that the time-derivative T a dtqr^ (of order e 2 ) as well as the 
nonlinear terms (of order e 2 at most) give negligible contributions. Up to 
the order e 2 , the final moment equations read [TJ] 

n a T a e rmn q^b n + r a G?W + + = o + 0(e 2 ) 
n a T a e rmn q^% n + Ta G°® + 5 ? (3) + 9? {3) = + 0(e 2 ) (3) 
n a T a e rmn q^b n + r a G a r ^ + = + 0(e 2 ) 

The moment equations link the generalized frictions Gr^ with the purely 
thermodynamic forces gr, g?^ and the generalized stresses g~r ^ 2n+l \ \, n 
is a unit vector along the magnetic field B i.e., b n = B n /B, e rmn is the 
completely antisymmetric Levi-Civita symbol and 0, a is the electron (a = e) 
or the ion (a = i) Larmor frequency. From Eq.Q we obtain 

N 

S e = qU 9 V + qf)gf) - qWfW + £ q fn + l)- g e(2n + l) 

n=l 

N 

= q f>)gf) + £ g <(2n+l)^(a«+l) (4) 
n=l 
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The expression for the entropy production is therefore an infinite bilinear 
form, truncated at n = N. We draw attention to the fact that expression Q 
is valid up to the leading order in [i = m e /rrii and, therefore, the ion-electron 
entropy production, giving contribution of relative order \i = m e /mj, has 
been neglected. This also explains why the expression for the electronic 
entropy production differs from the ionic one. In literature, expression Q is 
referred to as the quasi-thermodynamic form of the entropy production |12j . 
If we work in the twenty-one moment (21 M) approximation (i.e., N = 2), 
Eqs Ml) can be cast into the form 



For easy reference, we list here explicitly the relations between the dimen- 
sionless and the corresponding dimensional Hermitian moments: 



where m a and T a are respectively the mass and the temperature of species 
a and u denotes the centre- of -mass velocity of the plasma. Moreover, j r , 
and Q° indicate the electric current, the particle fluxes and the heat fluxes, 
respectively. is the dimensional fifth-order Hermitian moment corre- 
sponding to ^ ■ For completeness, we also report the relation between the 
pressure tensor tt^ s and the second-order tensor Hermitian moment qrs 
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The dimensionless generalized forces ^ 2n+1 ^ and the dimensionless gener- 



alized source terms are defined as 
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5 a Vr//,, 



2 r (Zf^ 1/2 
5 a \m, 



^^V s (n a Tj^))+V2T- 7 / 2 V s (Tl/ 2 C( 2 )) 



Indicating with e the absolute value of the charge of the electron and with Z 
the charge number of the ions, we have e a = — e for electrons and e a = +Ze 
for ions. E is the electric field of the plasma. In the local dynamical triad it 
is possible to show the validity of the following relations [12] : 



g^ n) = 0(e); = 0(e); #W = 0(e); g a A {n) = 0(e) 

L A = 0(e); L ± = 0(e 2 ) 



(9) 



where L is any dimensionless transport coefficient. Taking into account 
relation Q, expressions §5§ reduce to 



S e _ „(1) /•„(!) 



i(3)/ i(3) _i(3) 



_e(l)^ + „e(3)^e(.S> 



e(5) e(5) _ *e(l) (1) ~e(3) e ( 3 ) 



X* = gj w fo; w + gf) + 4 5) 4 5) + 4?sf > 



(10) 



Eqs {| 1 0|) are valid up to order e 2 . q e J^' , qffi and qty are the radial Her- 



mitian moments truncated at the order e 2 . Notice that = — q^ 
%d 9p^ an< ^ — Qpci9p are * ne classical contributions to the entropy 
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productions, can further be decomposed as 

(1) e(l)P e(l)$ e(l)A e(l) 

^ (1)P --f^) V2 -e— V„P e 

1 VT e / rn e n e 



e(l)* _ /"T-eX 1 / 2 e 



m e 

" V T e / m e II 

where P e is the electron pressure, — V r $ represents the self-consistent electric 
field built inside the plasma and Ej- represents the electric field induced 
by purely external means (such as transformer coils) in the confined plasma. 
The self-consistent electric field is, usually, much smaller than the external 
electric field and for this it can be neglected. 

It is useful to collect here the hypotheses adopted to obtain expressions 
(Hob: 



1) The state of the plasma is not too far from the reference local equilibrium 

state; 

2) The drift approximation is applicable; 

3) The ratio /x = m e /mj -C 1 and therefore quantities of order /j, have been 

neglected (but not quantities of order less that ji ); 

4) The tensor Hermitian moments have been neglected; 

5) The expressions have been evaluated in the 21 M approximation. 

Assumption 1) deserves an additional comment. The expression for the 
entropy production has been obtained making use of the so-called local equi- 
librium principle. However, the validity of this expression goes beyond the 



validity of the hypothesis of local equilibrium i.e., Eq. (10) remains correct 
also when the local equilibrium principle is invalid. Indeed, it is possible to 
show that the expression of the entropy production, written as sum of prod- 
ucts of thermodynamic forces and their conjugate fluxes, is valid throughout 
the whole range of thermodynamics [13]-[T3j. More generally, we can state 
that the limit of validity for the expression of the entropy production, writ- 
ten in a bilinear form, establishes the limit of validity of the thermodynamic 
description of a physical system [13] • In ref.[15] we can find many physical 
examples where the local equilibrium principle is violated and the expression 
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for the entropy production can still be brought in a bilinear form. Then, 
the only precaution that we have to take is to truncate the expressions at 
the first order of the drift parameter e. 

From the physical point of view, we have taken into account one assump- 
tion and the plasmadynamical balance equations: 

a) Collisions are the only source of irreversibility or dissipation; 

b) The plasmadynamical equations, expressing the conservation of mass, 

energy and momentum of plasmas, have been taken into account; 

c) The plasma is in mechanical equilibrium i.e., ^ = 0. 

Assumption a) is valid whenever the plasma is quiescent. Additional sources 
of dissipation must be taken into account if the plasma becomes unstable 
and turbulent [16]. The hydrodynamic equations have been used to obtain 
the Hermitian moment equations. As already mentioned, the truncation of 
the hydrodynamic equations is based on the smallness of the drift parameter 
e and not on the hydrodynamic parameter Ah, defined as the ratio between 
the shortest mean free path of the particles and the largest length of the 
hydrodynamic gradients. The momentum conservation implies the following 
relation 

ft(^) 3/ V + fi(f) 3 V» = (12) 
where Pi denotes the ion pressure. In terms of dimensionless generalized 



thermodynamic forces, Eq.(12) can be cast in the form |12j : 



5; (1) +a 5 ; (1) +^ (1) +^; (1) = o (is) 

wherea=^(^J . 

The Hermitian moments q» are linked to the poloidal fluxes cu® through 
the following relations |12j : 



u± = acu 1 — u>i 



(l) (i)p . B 



K a —gT> + —to? K a = (14) 
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where a = (fiTi/Te) 1 / 2 and Q a Q is the Larmor frequency associated with 

the average magnetic field. (3q =< B 2 >*/ 2 where < • • • > denotes the 

magnetic- surface averaging operation. B^ and Bg are the components of the 

magnetic field in the local geometric triad (often called the Hinton-Hazeltine 

coordinates [T7]) i.e., B = B$eg + B^e^. eg and are the unit vectors in 

(i)p 

the local geometrical triad. Quantity g p is defined as 

where P indicates the total pressure i.e., P = P e + P{. The poloidal fluxes 
are surface quantities, independent of the poloidal angle 6. In general, the 
quantity K a is not a surface quantity but it becomes one for some configu- 
rations of the magnetic field as, for example, in the standard model. In this 
paper we decide to make the assumption that also quantity K a is a surface 
quantity. 

The parallel fluxes q"^ and q^' can be further decomposed in three 
contributions 

a(n) a(n) , a(n) , a(n) 
9|| =%cl + 

„(1) _ (1) , >) , (1) (16) 
q \\ ~ q \\cl + q \\ps ^ q \\bp y L0 > 

The first contribution, <?juj , represents the parallel flux evaluated by the 

classical transport theory. The second contribution, q?J^ , is the parallel 
Pfirsch-Schluter fluxes, which we define as 

< Bqffi >= V a, n (17) 

The differences between the total parallel fluxes and the classical plus the 
Pfirsch-Schluter fluxes will be referred to as the banana/plateau flux q"f£ 



We can easily calculate the value w°* of the poloidal fluxes such that Eq.( 17) 



is verified. Indeed, from Eqs (14) and (16) we find 



<Bq^>=<B(-K e ^ p + ^u* 1 )>=0 <4 = K e gV> p 

< Bq«® >=< B {K a ^gf) + £o£*) >= ^ = -K^f* (18) 

< B(qf^) >=< b£u?) >= W5 a * = 



10 



From Eq.(14), we have 
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where 
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o(5) _ a(B) 
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' 1 1 6pcZ 
' 1 1 bpcl 



B 



a(n) 
1\\bp 

l\\bp 



'\\bpcl 



(19) 



+ 



+ 



a(n) 
'lie/ 



(20) 



and 



W = 

ps 

"(3) 
IIp s 

"(5) 
Hp* 



B__Po 
0o B 
B__fh 



K o (1)p 



K a 01 ^ 



(21) 







As we can see, using Eq.(17), we have easily obtained the expression of the 



parallel fluxes in the Pfirsch-Schluter regime. In literature, the contribution 
to electric current (3q < produced by the pressure and the thermal 

gradients, but not by the parallel electric field, is referred to as the bootstrap 
current [18], [19] and [20 . The presence of the bootstrap current leads to 
the attractive idea of operating a tokamak in a steady state, with E^ A ) = 0. 



From Eq.(17) we immediately found that there is no contribution to the 
bootstrap current in the Pfirsch-Schluter transport. On the other hand, 
the classical transport contributes to the parallel electric flux only through 
the presence of the parallel electric field. Therefore, the bootstrap current 



exists only in the banana regime. Taking into account Eqs (16) and (20) 



the expression of the entropy production, Eq.(10), can be brought into the 
following form 



-W („W _ n e ^>) + n e(0) (n e(0 > + rf w l 4- n [1 > (a [1> - a e{1 >) 

\\p S \9\\ 9\\ >^ q \\psyy\\ > + %b P cA9\\ 9\\ ) 



=e(lh 



e(3)/ e(3) _e(3) 



(!) iJX) 



e(3) , e(3) 

i\\bpd\9\ 



, ri.)j . -e(3)\ . 

(V +9\\ ) + Q 



e(5) -e(5) *e(l) n (l)p ^(3)^(3) 



\\bpd9\\ 



%ci9, 



+ 



1 P cl 9 P 



&^ + tf ) + fiCtf + ^f 3) ) + &S (5) + ^? 3) (22) 



l \\bpcl 



\\bpcl- 



where we have used the identity g p 



(i) 



-9p 1)p m 
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The quantities of major interest are the radial fluxes averaged over a 
magnetic surface < qp^ >. Indeed, these are the quantities measuring the 
leakage of matter and energy through the confinement region. It is possible 
to show that the average radial fluxes can be decomposed in six contributions 
[12] and [21] 

< qf n) >=< qf n) >d + < qf n) > P s + < qf n) >b + < qf n) > P 

+5 nA < > dr +5 nA < g"<"> > mdr (23) 

In general the electric drift fluxes < q p ^ >d- r and the modified electric drift 
fluxes S Uj i < qp > m dr are very small contributions compared to the other 
fluxes. Therefore, they will be neglected in the forthcoming sections. The 
total average radial fluxes are simply the sum of the classical, the Pfirsch- 
Schluter, the banana and the plateau averaged radial fluxes: 

< q^ >^< #W > d + < g;W > ps + < tfW > b +< > p (24) 

In section ^ we shall provide the nonlinear corrections to the averaged 
radial fluxes in the nonlinear classical and Pfirsch-Schluter regimes. 

Finally, it is important to note that our formalism starts from the concept 
of the entropy production. Thus when the plasma is far from equilibrium, 
the entropy production is no longer a potential and its extremal property 
does not determine the non- equilibrium steady states of the plasma. A cor- 
rect theory cannot be based on the minimum entropy production theorem 
for finding steady states when the plasma is far from equilibrium. In fact, 
a plasma far from equilibrium must satisfy the Universal Criterion of Evo- 
lution established in 1954 by Glansdorff and Prigogine [22]. This criterion 
however is not derived from a variational principle and by itself is not able 
to provide the corrections to the Onsager theory. The nonlinear corrections 
will be obtained by solving the field equations [1] combined with nonequilib- 
rium statistical mechanics, establishing in this way the link between micro- 
and macro-levels. In section Q we briefly summarize the Thermodynamic 
Field Theory and we write the field equations for a fully ionized plasma 
beyond the Onsager region. 

3 Linear Analysis 

In the linear region, the relations between the dimensionless Hermitian mo- 
ment and the dimensionless thermodynamic forces can be brought into the 
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following form: 

gJ (3) =a | |( 5 i 1) -^ (1) ) + K|(^+5; (3) ) (25) 
for the parallel electron fluxes, and 

^ = ^f l) +gf l) ) (26) 

for the parallel ion fluxes. Coefficients cry, ay, rejj are the dimensionless par- 
allel component of the electronic conductivity, the thermoelectric coefficient 
and the electric (i = e) or ion (i = i) thermal conductivity, respectively. 
The relations between the dimensional and the dimensionless transport co- 
efficients are 

_ e 2 n e _ 
m e 

On = W Te^n (27) 

V 2 m e 

q. 5 n a T a _ a 

2 m a 

The parallel transport coefficients define a definite positive matrix [12] 

or||, k|, k\ > 

or[|K| - a\ > (28) 



Eq.(25) must be supplemented by the solubility conditions expressing the 



mechanical equilibrium of the plasma: 

g^ + ag^+gf^agf^O (29) 

Up to the order e 2 , the classical electron and ion radial Hermitian moments 
are linked to the thermodynamic forces by the following relations 

C = <M 1)P - &±#*> + 0{e 2 ) 
C = -&±# )P + > + 0(e 2 ) (30) 

+ <>{?) 
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where cj_ , a± and kj_ are the dimensionless perpendicular component of the 
transport coefficients. These coefficients satisfy the relations 



a±, h e ± , Kj_ > 
a±K e ± — a\ > 



(31) 



The transport equations are derived by taking into account the definition 
of the fluxes in the different regimes. In the linear regime, the classical 
transport equations are immediately provided by Eq. ( 30 ) : 



< 



-e(l) 
Ipcl 



>=< <r± > g 



(i)p. 



< a± > 5p (3) 



< c > 



(32) 



< ^ >=< *i > a f ] 



Notice that the classical dimensionless radial fluxes can also be obtained by 
the following definition (see, for example refs [12 and [23J) 



< C ] > 



1 



a(n) 



> 



The dimensionless radial Pfirsch-Schluter fluxes are defined as 

i? 2 

PI 

From Eqs Q , ( 25 ) and ( 34 ) , we easily find 



< a(n) > _ 



-K a r a <^l\ 



Gf n) > 



< 



< 



< 



A 1 ) 



,i(3) 



MP 



<W 3) ) 



> ps =Ki{g- l)(c e n g p 
>ps =K 2 e (g-l)(-ci 3 gW p + 4 3 g; 



> 



K?(g-l)c 33 gl 



where factor g is defined as 



(33) 



(34) 



(35) 



g 



1% < & > 



and en, C13, C33 are the electron collision matrix elements 

K ll ce\\ Cm 



CllK' 



a; 



'13 



CllK 



II HI 



a: 



2 C 33 



Cm K„ 



a 



(36) 



(37) 



From now on, we shall adapt a more compact notation. We shall label 
with X% (fi = 1,2) the elect ron thermodynamic forces and with Y4 the 
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ion thermodynamic forces. Symbol A distinguishes the different regimes 
A = (cl, ps). The conjugate electron and ion thermodynamic fluxes will be 



denoted with symbols j\fl and J]l' , respectively. We have 



(0 



ps 



■I 



(e) 



fj,cl 



/(<0 




and the Onsager matrixes read 



L 



(e)clfiu 



L 



(e)pspv 





— a± 


\-a± 


hi 


f!" 


OL\\\ 







Y d 



Y, 



ps 



9 r 



i(3) . _i(3) 
9\\ +9\\ 



J. 



(<) 



5 i(3) 



V/3 



B 



L(i)d — ftj_ 



L{i)p S — ^ 1 1 



(38) 
(39) 
(40) 

^5? 3) (41) 

(42) 
(43) 



Our objective is to evaluate the nonlinear terms that should be added to 



the expressions (25), (26) and (30) when the plasma is far from equilibrium. 
These nonlinear corrections are provided by the solutions of the thermody- 
namic field equations. With these solutions, through the expressions (33) 



and (34), we shall provide the nonlinear transport equations in the differ- 



ent regimes. The definitions (33) and (34) for radial fluxes in the different 
regimes assumes ambipolar diffusion i.e., the electric current circulates only 
on the magnetic surfaces: < T e p >= Z < T p >. From the neoclassical the- 
ory we find that not only the total average fluxes, but also the separate 
classical, banana and Pfirsch-Schluter fluxes are ambipolar. The nonlinear 
corrections are determined by combining the solutions of the thermodynamic 



field equations with the expressions (33) and (34) and, therefore, the classi- 



cal and Pfirsch-Schluter fluxes (as well as the banana flux) remain ambipolar 
also in the nonlinear regime. 
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The Thermodynamic Field Equations for Mag- 
netically Confined Plasmas 



it is known that the validity of the Onsager reciprocity relations (25), (26) 



and (30) is limited to the so called linear region. As already mentioned in 
the introduction, a covariant thermodynamic field theory (TFT) has been 
proposed in 1999 in order to evaluate how the relations flux-forces deform 
when the thermodynamic system is far from the linear ("Onsager") region. 
Let us briefly summarize the main aspects of this theory when the transport 
coefficients are fully symmetric. 

Clearly, one cannot develop a transport theory without a knowledge of mi- 
croscopic dynamical laws. Transport theory is only but an aspect of non- 
equilibrium statistical mechanics, which provides the links between micro 
and macro-levels. This link appears indirectly in the "unperturbed" metric, 
i.e., the linear classical and neoclassical transport coefficients used as an 
input in the equations: these coefficients have to be calculated in the usual 
way by the methods of kinetic theory. In the TFT description, a thermody- 
namic configuration corresponds to a point in the thermodynamic space Ts, 
defined as an n-dimensional manifold covered by the n independent thermo- 
dynamic forces X^. The elements of the metric tensor g^ are the symmetric 
components of the transport coefficients t^ u . The conjugate flows, Ju, are 
dual to the thermodynamic forces through the relation = T ilv X v . In this 
equation, as in the remainder of this paper, the Einstein summation con- 
vention on the repeated indexes is adopted. The positive defmiteness of the 
matrix g^ v ensures the validity of the second principle of thermodynamics 
[1] . The theory is constructed with the following requirements ( [T] , [26] ) : 

• The theorems valid when a generic thermodynamic system is out from 
equilibrium should be satisfied] 



The thermodynamic field equations should be invariant under the ther- 
modynamic coordinate transformations (TCT). 



With the elements of the transport coefficients we construct two objects: 
operators, which may act on thermodynamic tensorial objects and thermo- 
dynamic tensorial objects, which under coordinate (forces) transformations, 
obey well specified transformation rules. 

Operators 

We introduce two operators, the entropy production operator a and the 
dissipative quantity operator P{X) considering the transport coefficients as 



16 



elements of n x n matrixes r and g, which multiply the thermodynamic 
forces X expressed as n x 1 column matrices: 



a(X) a(X) = XgX T 

P(X) :-► P(X) = XtX t (44) 
[The dot symbol stands for derivative with respect to parameter 77, defined 



in Eq. (49)]. Thermodynamic states X s such that 

P{X a ) = (45) 
are referred to as steady- states. These are physical quantities and should re- 



main invariant under thermodynamic coordinate transformations. Eqs (44) 
should not be interpreted as the metric tensor g^ u , which acts on the coor- 
dinates. The metric tensor acts only on elements of the tangent space (like 
dX^) or on the thermodynamic tensorial objects. 

Transformation Rules of Entropy Production, Forces, Flows and 
Thermodynamic Tensorial Objects 

According to Prigogine's statement [13 , thermodynamic systems are 
thermodynamically equivalent if, under transformation of fluxes and forces, 
the bilinear form of the entropy production a remains unaltered. In mathe- 
matical terms, this implies: 

a = = J'^X'» (46) 

This condition requires that the transformed thermodynamic forces and 
flows satisfy the relation 

These transformations are referred to as Thermodynamic Coordinate Trans- 
formations (TCT). 

By direct inspection, it is easy to verify that the general solutions of equa- 



tions (47) are 



X '» = X 1 F»(— — ■■■ XU \ (48) 

where is an arbitrary function of variables X^ jX^ with (j = 2, . . . , n). 
We immediately note that dX^ and d/dX^ transform like a contra- variant 
and a covariant vector, respectively. Vectors dX^ define then the tangent 
space to Ts. It also follows that the operator P(X), i.e. the dissipation 
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quantity, and in particular the definition of steady-states, are invariant under 
TCT. Moreover, it is easy to check that the metric tensor g^ v transforms like 
a thermodynamic tensor of second rank. In particular, parameter rj, defined 
as 

dr, 2 = g^dX»dX v (49) 

is a scalar under TCT. 
Operator O 

is also invariant under TCT. This operator plays an important role in the 
formalism. 

The Principle of Least Action 

We introduce now the following postulate: 
There exists a thermodynamic action I, scalar under TCT, which is sta- 
tionary with respect to arbitrary variations in g^ v . 

We want to construct an action, scalar under TCT, from the metric tensor 
and its first and second derivatives, that has linear second derivatives. It is 
easily checked that the only action with these requirements is ([I] and |26j ) 

/ = J (iv<T + KX^x" ^ d Qxv ) dV = J(R + ^o(gx K )0(g XK ))dv 

(51) 

where 

A (j = constant having dimension the inverse of the entropy production 
R»u = rj v - rJ„ iA + r^r^ - r^rJ A (The Ricci Tensor) 
T„ v = affine connection 

dV = infinitesimal volume element of Ts (52) 



To avoid misunderstanding, while it is correct to mention that this postulate 
affirms the possibility of deriving the thermodynamic field equations by a 
variational principle it does not state that the solution of the thermodynamic 
field equations can also be derived by a variational principle. In particular 
the well-known Universal Criterion of Evolution established by Glansdorff- 
Prigogine cannot be derived by a variational principle. 
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The Thermodynamic Field Equations 

We immagine g^ u to be subject to an infinitesimal variation g^ u — > g^ u + 
fig^u where 5g^ u is arbitrary, except that it is required to vanish as | X^ \— > 
oc. By imposing that the action (51 ) is stationary with respect to arbitrary 
variations in g^, we find 

+ \ a S» v = (53) 

where 

= -20\gn+^0{gnO{g uX ) + l -g^O{g XK )0{g^) 
-[0{\ng)+A]0{gn ' (54) 

Close to the Onsager region, having metric L^v, we can write 
9iiu = L^ u + h^v + 0(e 2 ) 

X a = 0{e) with e = Ma X { El9enValue f 9 ^- L ^ ] \<<l( 55 ) 

I Eigenvalues[L^\ ) 

and therefore, h^ u are small variations with respect to Onsager's coefficients 
satisfying the equations [1^ 

t Xk d 2 h^ Xk d 2 h XK Xk d 2 h Xu Xk d 2 h x ^ = 

dX x dX K dX^dX" dX K dX^ dX K dX» { ' 



Eqs (56) should be solved with the appropriate gauge-choice and boundary 
conditions. 

Looking at the expression of the entropy production for electrons and ions, 
Eq.(lO), we immediately realize that the thermodynamic space of a fully 
ionized plasma possesses 11-dimensions. Notice that in this case, the coef- 
ficients of the reciprocity relations are fully symmetric and, therefore, the 
metric tensor coefficients g^ v coincide with the transport coefficients of the 
plasma. Luckily, the thermodynamic space results to be the direct sum 
of six independent subspaces: three subspaces for electrons and three sub- 
spaces for ions. The electron subspaces have dimensions 2 for the classical 
and Pfirsch-Schliiter transport and 3 for the plateau/banana regimes. The 
ionic subspaces have dimensions 1 for the classical and the Pfirsch-Schliiter 
transport and 2 for the banana regime. Below, we report the schematic 
decomposition of the thermodynamic space: 

Ts 11 = Ts\ e) © Tsfa where 

Ts \e) = Ts d(e) © Ts ps(e) © Ts b/p{e) 

Ts \i) = Ts \i) © Ts )>s(i) © Ts b/p(i) ( 57 ) 
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In sections ^ and ^ we have found the expressions for the entropy of a 
fully ionized plasma and the linear relations between the thermodynamic 
forces and the conjugate fluxes. 



Considering Eq.(lO), the electron and ion weak- field equations read 



ffl h (e)A & h (e)A &h& A d 2 h {e)A 

jXk d h »V , jXk d h \K jXk g h \v J-XK Xfl » 

^ A dx\ex\ ^ A dx\dx\ ^ A dx\dx\ ^ A dx\dx\ 
f)2 h (i)A & h m ^ h m a 2 h { ! )A 

j Xk d h VV , j Xk 6 h X K __ r Xk 6 h Xu _ r Xk _ Xfi 

® A dY*dY* ® A dY£dYX ^ A dYldY» ^ A dY£8YX 

(58) 

where A = (cl, ps, p, b), corresponds to classical (cl), Pfirsch-Schliiter 



(ps), plateau (p) and banana (b) regimes. To solve Eqs (58), we have to 
chose the particular system where the calculations [i.e., the gauge choice, 
see appendix will be performed and then to solve the resulting P.D.E. 
with the appropriate boundary conditions. These tasks will be accomplished 
in the next sections. Notice that, it follows from the second principle of 
thermodynamics that these P.D.E. are of elliptic type. 



The Electron and Ion Gauge-invariant Solutions 
for the Collisional Regimes 



To solve Eqs (58) it is convenient to work in the harmonic coordinates. We 
label X^ (/i = 1,2 and A = cl,ps) the electron thermodynamic forces and 
with hffi d (X c i) and hffi ps (X ps ) the small field perturbations h^ u when we 
refer to the classical and Pfirsch-Schliiter electron regimes, respectively. In 
harmonic coordinates these fields satisfy the following equations 

Xv _ Ij-^A a V 
or(e)ps -, f)h^ ps 

r^x ori Xv _ 1 A Va , . 

(e)P s dX» s 2 dX U s W 

where and indicate the corresponding electron Onsagers's ma- 

trices in the classical and Pfirsch-Schliiter regimes, respectively. In the 
next section is shown that in the weak field approximation, the ion clas- 
sical and Pfirsch-Schliiter regimes remain unperturbed (i.e., hfy A (YA) = 0). 
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If h$ A (X A ) are cast into the form h$ A {X A ) = c$ A f^ A (X A ) with c$ A 
constant and f^ A (X A ) functions depending on the electron thermodynamic 
forces, Eqs ( [59] ) give 

(e)d 1 T \ K (e)d,A df(X d ) _ 

L {e)cl C \v ~ 2 (e)ci C A K °v J ~QX^~ ~ 

\ df(x ps ) _ 



J {e)ps L '\v 2 ( e )P s Xk v ) dX 

/ t 

which are verified for any value of ^^^^ om y if 



ps 



C 



( e ) c ' — -v T 
C\nF = XpsL^pg^u (61) 

whit Xd and Xps being parameters independent of the thermodynamic forces. 

6 The Nonlinear Classical and Pfirsch - Schliiter 
Transport Regimes 

In the nonlinear region, the relations between the dimensionless Hermitian 
moments and the dimensionless thermodynamic forces can be cast into the 
form: 

& = (1 + Xdf {e)d (X d )) ( a ±9 U p - a^gf) 



pel 

= (! + Xd& )d {X d )) {-a x9 ? P + «iflf>J (62) 

C = («i + 4f )4 3) 

for the classical regime, and 

?S 1J = (1 +Xpsf^ ps (X ps ))(a {l (g^ -g^) + a^gf } +g°®j) 

if* = (1 + Xpsf^iX^a^ - gf>) + R\{g^ + gf)) (63) 

?i (3) =(4+4? s )(5j (3) +5j (3) ) 



for the Pfirsch - Schliiter regime. In these equations, we have taken into 
account that the perturbations h$, A are of the form 

h$ A {X A ) = XAL {e)A ^f {e)A (X A ) (with A = cl,ps) (64) 
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, e) by solving the field equations. In 



Let us now determine h^J (a 
harmonic coordinates, the ion field equations for the classical and Pfirsch- 
Schliiter regimes take the form 



r-l 

^(i)cl 



d 2 h 



(i)cl 
33 



8*1 

2ui i )P s 







L 



d 2 h 



33 



(i)ps Qy2 
w 1 ps 



(65) 



The solutions of Eqs (65) are 



h 



(i)cl . 
l 33 " 
(i)ps 
33 



-- a 2 Yp S + 6 2 



(66) 



where a n and b n are constant coefficients. The nonlinear solutions must co- 
incide with the Onsager solutions when the system approaches equilibrium. 
This condition imposes that b\ = b 2 = 0. Moreover, solutions (66) must 
verify the harmonic gauge conditions: 



L 



-1 



dh\ 



(i)cl 



33 



1 



L 



-1 



dh 



(i)cl 
33 



L 



dh 



(i)ps 
33 



(i)ps qy 



ps 



I L -1 ^33 

2 Mp s ay r 



(i)ps 



(67) 



Inserting solutions (66) into Eqs (67), we find a\ = a 2 = 0. Therefore, the 



Thermodynamic Field Theory does not correct the expressions provided by 
the neoclassical transport theory in the classical and Pfirsch-Schliiter trans- 
port. Notice that we could have obtained the same results simply by observ- 
ing that in the classical and Pfirsch-Schliiter regimes the ion field equations 



[i.e., the second equation in Eqs (58)] vanish identically. 
The electric field equations are P.D.E. of elliptic type, 



dh 



(e)c2 

flU 



(e)dQ X »Q X v 







dh 



(e)ps 

fJ,U 



(e)psQ X g s Q X V s 







(68) 



These equations can be reduced to Laplacian equations 



= 



(69) 



22 



by performing the following coordinate transformations 



x d 
xl 

v-l 
^ps 
v2 
-^ps 



x\ d 

X2 d 

Xlpg 

X2p S 



(70) 



The "scaled" variables XI, 2 d and Xl,2 ps are defined as follows 



^ + cl 

X2 cl 

> X -d / 
' Xlp s ' 

ps 
'ps 
ps , 



u 



x d 
xl 



A 

x% 

.A 



u 



■^ps 
ps 1 y2 



(71) 



where \± d and U c / are, respectively, the eigenvalues and the matrix formed 
by the eigenvectors u +d ,U- d (put in column) of the contravariant Onsager 
matrix L^ c[ i.e., 



^(e)cl 



1 



K c ± a± — a\ \a± a± 



a± 



with A 



2{R e L a± — a\) 
U d = (u +d ,u^ d ) where 

V2a ± 



d 



2a ± 
1 



(72) 



The corresponding relations for the Pfirsch-Schluter regime are 
1 



[e)ps 



A± 



<7|| K| — CK'| 

0T|| + fefj ± yA 



— a | 

CT|| 



ps 



ps 



2(a»K 



with A 



ps 



U ps = (u H 



where 



^) 2 + 4a 2 



cr II — sJ±i/A t 



(73) 



ay - Sp 2 + 4a 2 ± (ay - fcjj) y/A 



■ps 



2a„ 
1 



The solutions of Eq.(69) are submitted to the following conditions: 
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a) The solutions vanish at the origin of the axes ; 



b) The solutions are constant, with value Xci^(e)d^v m t ne classical regime, 

and XpsL[ e ) ps ^ v in the Pfirsch-Schliiter regime, on the circumferences 
with radii equal to Rod and Ro ps , respectively; 

c) The solutions are invariant with respect to the permutation of the axes 

Xl c i and X2 c i for the first equation (69) and of the axes Xl ps and 
X2 ps for the second equation (69); 



d) Inside the circumferences, the solutions are of class C 2 (possibly, except 

in a set of measure zero where they should be at least of class C°) 

e) The profiles for the losses (matter and energy) are, at least, functions of 

class C°. 

The first condition is due to the fact that our solution must coincide with the 
Onsager matrix when the system approaches equilibrium. Point b) provides 



the boundary condition of Eqs (69) at thelimit of a bounded region. The 



existence of a bounded region, where we have to solve Eqs (69 ), is warranted 
by the fact that the amplitudes of the thermodynamic forces are, in 
reality, bounded (see figs [6jf|8] for a JET plasma). Therefore, 3Ma G 1Z + 
Ma = max{X^, fj, = 1,2} < oo. Since the thermodynamic forces are 
bounded and independent, symmetry implies that they span a (finite) ellipse 
in the thermodynamic space. On the other hand, it is always possible to 



perform a coordinate transformation ( 70 ) such that the thermodynamic field 



equations take the symmetric form of Eqs (69) and the ellipse transforms 
onto a circle of radius Rq. After the coordinate transformation (70 ), the new 
thermodynamic forces become perfectly equivalent while at the boundary, 
one direction cannot be privileged with respect to another. This condition 
requires that the solution is constant at the boundary. The third condition 
ensures that the two axes are treated in an equivalent manner. The next to 
last condition requires that the only solutions, acceptable from the physical 
point of view, are functions continuous in all points inside the circle (i.e., 



excluding the boundary) and, which are solutions of Eqs (69) (possibly, 
except in a set of measure zero). Finally, condition e) requires that the 
losses profiles are continuous functions. 

The first task is to construct a well-posed Dirichlet problem from the above- 
mentioned conditions. This has been accomplished in appendix In 
appendix ^ it is shown that, up to the first order, i.e., in the limit of 
validity of the weak-field approximation, the boundary conditions derived 
by conditions a)-e), respect the principle of covariance. Successively, we have 
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to find the solutions of the P.D.E. (69) submitted to the Dirichlet problem. 
These are obtained in appendix ( |9.1[ ) as 



h$ cl (Xl d ,X2 cl ) = XciL {e)dflu f^ cl (Xl d ,X2 cl ) 
hfy ps (Xl ps ,X2 ps ) = -L^ psflu f^ ps (Xl ps , X2 P 



(74) 
(75) 



with 



f^ d (Xl d ,X2 d 
f^ ps (Xl ps ,X2 



arctan 



7T 



Xl d X2 



cl 



Rt l -^il + x2lY 



arctan 



7T 



L cl ^ ^"d> 
Xlr,f,X2 r 



K>4 

-"■Ops 



(Xll + X2l s f 



(76) 



Notice that, \ d is still undetermined at this stage while the general condi- 
tion e) imposes Xps = — 1. Parameter Xd is easily determined by imposing 
that, for large values of the thermodynamic forces, the expressions for the 
diffusion coefficient evaluated by the TFT and the kinetic theory must co- 
incide. Through a simple kinetic model, it is possible to show that Xd = 1 

mi- 

In the classical regime, the nonlinear transport equations are immediately 
obtained from Eqs ( 62 ) . They read 

~e(l) ( , ~ 2 

>= ( (<<7j_> + — <a± arctan 

7T 



< 



'pel 



Xl d X2 



cl 



<a± >+— <a± arctan 

7T 



R ii- 



(Xl 2 d + X2lY 
I Xl d X2 d | 



(i)p 



R 0cl 



{Xll + X2lf 



< 



Ipcl 



<a± >+— <aj_ arctan 

7T 



>)rf 3) 



Rt c i-^il + X2lY 



< 



+ ( <«I >H — < k± arctan 

V 7T 

5 *(3) ~i i( 3 ) 



L cl 

X\ d X2 d 



>)# )P 



(X11 + X2W 



>)9f ) 



(77) 



where XI, 2 d are given by Eqs (71 ) and (72 ). The nonlinear Pfosch-Schluter 
transport equations are obtained recalling Eqs ([3| and (34). They suggest 
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that 



^ Hp ^ps~ 



K R < 



B 



< 



A 3 ) 



> ps =K e <^ll 



_ & 

PI 
B 2 



Si) 



9f ] ) > 



ps~ 



1 



P 2 o 
Pi 



+^ (3) )> 



(78) 



which, combined with Eqs (63), give 



< of) > ps = Ki{~g ps - 1) U l9 W p - cfo 



_e(3) 



< <?f > > ps = K £ 2 (5 PS - 1) (-c e 13 gV p + 4s9f } 



(79) 



A 3 ) 



<q^> > ps =Kf(g-l)4 s gf^ 

where g and the collision matrix elements Cp V are given by Eq. (36) and 
Eq. (37) respectively and g ps is defined as 

>-f)(IM 



i+ < 



arctan 



ARl vs \Xl ps X2 va \ 
Bi ps ~(Xll s +X2l s y 



(80) 



with Xl ps and X2 ps given by Eqs (71) and (73). 
It is useful to re- write Eq. ( 79 ) as follows 



< <£« > ps = T vs Kl{g - 1) (c\ ig W P - c\ 3 gf>) 



< 



< 



/(3) 
A 3 ) 



> ps = ^ s i^( 5 - 1) ( -c^W p + c^f) 



(81) 



> 



Kf(g-l)cl 3 g 



LA 3 ) 



where we have introduced the amplification factor T vs defined as 

9 P s - 1 



ps 



5-1 



(82) 



Note that the amplification factor T vs is equal to 1 when the perturbation 
function f^ s (X ps ) is set to zero and it may assume very large values as 



P 2 

-"-Ops- 
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An estimation of the radii Rq c i and Rq ps can be provided. We define with Lh 
the minimum length of the gradient of the macroscopic quantities, namely 
density, the velocity, the pressure or the temperature: 



1 AJ Vi I 

— = Max '. v / 
L H A(X) 



(83) 



Without going into more accurate estimates, we note that these gradients 
pertain to macroscopic devices. Thus, the length of the gradients is very 
large compared to molecular dimensions. From Eqs ([8| and ( 15 ) we find 
after a little algebra: 



Ro d = Max{(Xl d + X2 d ) 1/2 } 



(84) 



Max< 



-{(-, 



T, \ 1/2 



a 1 + 



Pp 



+ V10 6 1 + 



Pi, 
P P 



5 ■ 
+ 2 C 



1/2' 



where 



a = a+Z? 2 . + a-/? 2 . 
1 



with 



a± 



Max 



P 



b = a + f3 + + a_/3_ 

} 



a+ + q;_ 



V p T e 



(85) 



o 



2 [^ ± 2 + o\ + 2ai ± («i + a ± )v^] 



R e L -d±± \fKa 



2a 



A, 



a±) 2 + 4ai 



Therefore we find 



TFT 



Neocl. 



(86) 



(87) 



i.e., the value of Rq c i evaluated by the TFT coincides with that one estimated 
by the neoclassical theory. Similarly, from Eqs Q, (15) and (21) we obtain, 
for Z = 1, 

R 0ps = Max{(Xl ps + X2 ps ) 1 ^ 
00 



-^-Max<JVp S | K e | T e — 



B 2 



0l 



Te. \ V2 T 

m, 



0.932 1 + 



+ 2.331 



1/2' 



The relation between the TFT and the neoclassical theory estimations is 
then 

Pops \TFT = ^°P S I Neocl 

Max{^ ps } (89) 
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The exact values of the parameters RqcL and Ro ps can be obtained numer- 
ically if the magnetic configuration and the profiles of the thermodynamic 
forces are specified. 

Fig.Q depicts, in polar coordinates, the function f^ A [XjC) (with A = 
cl,ps). We note that, in the weak-field approximation, the TFT does not 
correct the pure effects i.e., it does not correct the pure Fick's law and the 
pure Fourier's law. The TFT provides nonlinear corrections only for the 
cross-effects. At the end of this section we shall see that, for a JET plasma 




Figure 1: /^ a (Xa) (with A = cl,ps) in polar coordinates. This perturba- 
tion function satisfies the condition < f^ e ' A (Xj\) < 1. 

in L-mode, the fluxes of electron energy (heat flow) is further enhanced 
in the nonlinear classical and Pfirsch-Schluter (P-S) transports. The main 
conclusion of the present analysis is thus: 

1) In the nonlinear regime, the transport coefficients are given by the neo- 
classical values supplemented by extra terms evaluated by the TFT. 
These corrective terms are given by an amplification function times the 
Onsager matrix. This function corresponds to the arc-tangent function 
of the thermodynamic forces for the classical regime and to the ampli- 
fication factor T vs for the Pfirsch-Schluter transport. We have already 
mentioned that, through a kinetic model introduced in ref. it is 

possible to show that, in the collisional transport regimes, the electric 
diffusion coefficient and the electrical thermal coefficient are approx- 
imatively given by the linear (Onsager) transport coefficients times 



2S 



a universal function expressed by the inverse of the electron collision 
time r e . Since, the collision time decreases as the intensity of the ther- 
modynamic forces increases, the model predicts that the values of the 
electrical thermal and diffusion coefficients increase as the intensity of 
the thermodynamic forces increases; 

2) In the classical and Pfirsch-Schliiter regimes, the TFT does not correct 

the expressions of the ion heat fluxes evaluated by the neoclassical 
theory; 

3) The electron nonlinear classical transport coefficients may exceed the 

linear ones by a factor, which may be of order \d + 1 = 2. The 
electron nonlinear Pfirsch-Schliiter transport coefficients exceed the 
values evaluated by the neoclassical theory by the factor J- ps . Concrete 
calculations, for a JET plasma in L-mode, show that the amplification 
factor may assume values of order 10 2 (see below). These nonlinear 
corrections refer to the cross-effects; the pure Fick's and Fourier's laws 
remain, however, unaltered. 

This dissymmetry between the P-S ion and electron transport coefficients is 
encouraging since it is in qualitative agreement with experiment. Fig. ^ 
schematically illustrates the physical interpretation of the nonlinear Pfirsch- 
Schliiter effect [TT]. As the thermodynamic forces increase their magnitudes, 




Figure 2: The nonlinear Pfirsch-Schliiter effect. As the thermodynamic 
forces increase their magnitudes along the major radius, the resulting 
Pfirsch-Schliiter current must also increase. The intensity of this current 
is maximized when the thermodynamic forces reach their maxima values 
(i.e., Xl^ s + A2 2 s — > Rq ps ). 

after some collisions, the plasma pressure, through toroidal geometry, tends 
to be larger at larger major radii. The resulting increased pressure must be 
balanced by the internal magnetic force. As a consequence, the currents, 
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which produce this force, must increase their magnitudes producing an in- 
crement of charge accumulation in the upper and lower parts of the plasma. 
Charge accumulation is prevented by an increment of the intensity of the 
Pfirsch-Schluter current. Consistent with intuitive expectations, this effect 
is maximized when the thermodynamic forces reach their maxima values 
(Xl^ s + X2 2 vs 

loss (heat flux) and matter (particles flux) is possible. 
We conclude this section reporting the comparison between the neoclassical 
results and the TFT predictions for the electron mass flux and electron heat 
flow for L-mode JET plasmas. The numerical estimations of the electron 



Rqps)- In this limiting scenario, a large amount of energy 



fluxes, given by Eqs (78) and (79), can be evaluated by specifying the mag- 



netic field configuration and the expressions for the thermodynamic forces. 
In the local triad (e r ,Bg,B^), the magnetic field, in the standard hight as- 
pect ratio, low (3, circular tokamak equilibrium model, here referred to as 
the standard model, reads 



B = B G(r) b 9 + 



Da 



1 + (r/R) cos 9 



3 C 



with G[f) 



Rq(r) 



(90) 



where -Bo is a constant having the dimension of a magnetic field, q(r) denotes 
the safety factor and R is the major radius. For JET, Bn = 34500 Gauss and 
the major and the minor radii of the reactor are, respectively, R = 296 cm 
and a = 125 cm. The experimental profile of the q(r), for a JET plasma in 
L-mode, can be found in fig. [3] The specification of thermodynamic forces 
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Figure 3: The experimental safety factor versus the re-normalized minor 
radius p for L-mode JET plasmas [27] . 



g^ P and g"^ requires the knowledge of the electron and ion temperature 
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and density profiles. Figs [4] and [5] show typical experimental profiles of the 
electron temperature and the electron density of a L-mode JET plasms[^27 . 
In the standard configuration, the averaging formula of a physical quantity 




Figure 4: Typical experimental electron temperature profile versus the re- 
normalized minor radius p for L-mode JET plasmas |27] . 




Figure 5: Typical experimental electron density profile versus the re- 
normalize minor radius p for L-mode JET plasmas |27j . 

A(r, 9) simply reads [12] 

<A{r,9)>=^-[ d6(l + cosO) A(r, 6) (91) 
2-ir J \ R J 

Figures [6] [7] and [8] report the profiles of the thermodynamic forces for a 
L-mode JET plasma against the re-normalized minor radius p = -. In these 

2 At this stage, we can assume that the ion temperature profile coincides with the 
electron one. Notice that, for Z — 1, the electroneutrality condition imposes that ri ~ n e . 
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calculations Z = 1. The theoretical predictions of the Thermodynamic 
Field Theory (TFT) and the neoclassical results, for a JET plasma in the 
two collisional transport regimes, are compared in figs [9[rjl2| In particu- 



lar, Figs M and 10 report the (radial) electron fluxes evaluated by the TFT 



and the neoclassical theory, for plasmas in the classical and Pfirsch-Schluter 
regimes, respectively. The profiles of the (radial) electron heat fluxes, de- 
termined by these two theories in the classical and Pfirsch-Schluter regimes, 



are illustrated in figs 11 and 12 respectively. It is seen that, in the core of 



the plasma, the nonlinear classical transport coefficients exceed the linear 
ones by a factor, which may be of order 2. The electron nonlinear Pfirsch- 
Schluter transport coefficients exceed the linear ones by a factor, which may 
be of order 10 2 . Although having similar shapes, the curves strongly differ 



in magnitude (see Figs 14 and 15). 



(x 10 b ) 




(1)P 

Figure 6: The dimensionless thermodynamic force g p versus the re- 
normalized minor radius p for a L-mode JET plasma. This force is pro- 
portional to the gradient of the total pressure. 



7 Conclusions 

We apply the thermodynamic field to study the classical and Pfirsch-Schluter 
(P-S) transport regimes in a magnetically confined plasma. A new set of 
nonlinear transport equations i.e., the flux-force relations, have been de- 
rived. These equations provide nonlinear corrections to the linear ("On- 
sager") transport coefficients. A quite encouraging result is the fact that 
a dissymmetry appears between the P-S ion and electron transport coeffi- 
cients: the latter submits to a nonlinear correction, which makes the radial 
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g^ (3) (x io 6 ) 




Figure 7: The dimensionless thermodynamic force g p versus the re- 
normalized minor radius p for a L-mode JET plasma. This force is pro- 
portional to the electron temperature gradient. 



electron coefficients much larger than the former. Such a correction is absent 
for ions. This is in qualitative agreement with experiments. In particular, we 
have shown that when a plasma is out of the linear region, the classical and 
the Pfirsch-Schliiter transport coefficients are corrected by a universal func- 
tion times the Onsager matrix. For the classical regime, the amplification 
function is (approximatively) expressed as the arc-tangent function of the 
thermodynamic forces while in case of Pfirsch-Schliiter regime it corresponds 
to the amplification factor T ps . The theoretical predictions of TFT can be in- 
terpreted in terms of collisional mechanisms through heuristic kinetic model 
Through such a model, a crude estimation of the expressions for the 
transport coefficients can be derived. We showed that these coefficients are 
approximatively given by the linear (Onsager) transport coefficients times 
an amplification factor expressed as 

F P s = °4^^ (92) 
a eff. Te P s 

where T eps and 0" e //. are, respectively, the electron collision time and the 
effective electron cross section in the P-S regime and r e o and o e $fa are the 
electron collision time and the effective electron cross section estimated by 
the linear theory in this regime. In other words, the extra term found by 
the TFT is linked to an estimation of the time-step of the random walk of 
collisional particles. Combining our results with the kinetic model we find 

g-l /<T e //.o\ 2 



Te P s^^—( -^) r e0 (93) 
g ps — L \ a e ff. 
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Figure 8: The dimensionless thermodynamic force g p versus the re- 
normalized minor radius p for a L-mode JET plasma. This force is pro- 
portional to the ion temperature gradient. 

The electronic collision time decreases as the thermodynamic forces increase 
in magnitude. This implies that the electrical conductivity decreases as the 
intensity of the electric field increases, while the values of the electrical ther- 
mal and diffusion coefficients increase as the intensity of the thermodynamic 
forces increases jllj . These results are in line with numerical simulations of 
a Lorentz gas subject to a strong electric field [28] . 

Concrete calculations for L-mode JET plasmas have also been reported. 
We found that, in the core of the plasma, the nonlinear classical transport 
coefficients may exceed the linear classical ones by a factor of order 2. In 
this region of the plasma, for Z = 1, the electron nonlinear Pfirsch-Schluter 
coefficients may exceed the linear ones by a factor of order 10 2 . 
This analysis is extended to the banana and plateau regimes in part II of 
this article. 



8 Appendix: The Gauge Invariance of the Field 
Equations 

In the weak-field approximation, the small corrections h^ u to the Onsager 
matrix, i.e., g^ u (X) = Lp V + hfj_ v (X), satisfy the equations: 

j \k d 2 /y Xk d 2 h\ K Xk d 2 h Xu Xk d 2 h XfJL 

ax*x« ^ dxnx» dx^xv dx*x» { ' 

These equations should be solved with the appropriate boundary conditions. 
These will be explicitly derived in appendix Q . In this section we clarify two 
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Figure 9: L-mode JET plasmas, with Z = 1, in the classical regime. Predic- 
tions of the radial electron fluxes profiles according to the Thermodynamic 
Field Theory (TFT) and the neoclassical theory. 



important aspects: we provide the definition of gauge invariance of the field 
equations related to our problem and we prove that our choice of boundary 
conditions respects the principle of covariance. 

The most general coordinate transformation that leaves the field weak is of 
the form 

X » _> X'v = X^ + <?{X) (95) 

where ^g^P is at most of the same order of magnitude as h^, v {X). From 
Eq.(95), we have 



dX" 
dX u 

dX'v 



8£ + 



de^(X) 

dX» 
de v {X) 

dX^ 



(96) 



we also find 



9'^{X) = gU x ') + 



d9 §0 1 e\x) + o(<*) 



(97) 



and from the relation 



gU x> ) 



dX x dX K 
dX'v dX w 



gUx) 



(98) 
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Figure 10: L-mode JET plasmas, with Z = 1, in the Pfirsch-Schliiter 
regime. Predictions of the radial electron fluxes profiles according to the 
Thermodynamic Field Theory (TFT) and the neoclassical theory. 



we obtain 



9uu( x ') = 9tiv{X) - g\ u {X) 



de x (X) 



axM) 3 ^ 1 + 0{e 2 ) (99) 



Taking into account Eq. ( 98 ) , Eq. ( 97 ) becomes 



9L(X) = g^{X)-g\ u {X) 



From Eq.(lOO) we find 



h'{X) 



de x {X) 
dX^ 



, vx de x {X) , 

-gx^x) gj{u +- 



' 9 -^e\X) (100) 



de^X) de^jX) 
dX v dXv 



(101) 



The meaning of Eq. (101) is that the effect of an infinitesimal coordinate 



transformation on the tensor h^ v {X) is that the new tensor h'^^X) equals 

the tensor h fMU (X) — 9e Q X ^ ~ ^ax^ a * ^ e same coordinate point. 

Let us now consider the boundary conditions. As already stressed, the 
amplitudes of the thermodynamic forces are, in reality, bounded (see for 



example, figs Knf|8^ for a JET plasma). Therefore Eq.(94) should be solved 



with the appropriated boundary conditions. The choice of the boundary 
conditions cannot be arbitrary but should be done in order to respect the 
principle of covariance (other than the symmetry of the problem). In this 
paper we shall consider the case where, in coordinate X^, the boundary d 
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Figure 11: L-mode JET plasmas, with Z = 1, in the classical regime. 
Predictions of the radial electron heat fluxes profiles according to the TFT 
and the neoclassical theory. 



is an hyper-sphere B(Rq) of radius Rq, centered in the origin of the axes. 
h fll/ (X) equals a constant, say with value c^, on this hyper-sphere. This 
case corresponds to the situation in which in the thermodynamic space one 
direction is not privileged with respect to another one. In order to re-obtain 
the Onsager laws when the system approaches equilibrium, we have also to 
require that the tensor vanishes at the origin of the axes i.e., /tm,(0) = 0. 
In appendix Q, we show that these constrains, together with two other 
conditions, specify a well-posed Dirichlet's problem. 

We observe that, if h^ u (X) satisfies Eq.(94) and the boundary conditions, 
so will also be h'(X) = h tlv {X) — 9e gx^ — W5tP » obtained by performing a 

9e "^-%#^ vanishing 



general coordinate transformation, with the tensor 
on the hyper-sphere and at the origin of the axes. 
In conclusion, The coordinate transformation 



dx v 



X » _> X'" = X fl + e»(X) 



(102) 



with 



de^X) de u (X) 



dX v 



dX^ 



de^(X) de v {X) 



dX v 



dX^ 



B{R ) 



X=0 



(103) 
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Figure 12: L-mode JET plasmas, with Z = 1, in the Pfirsch-Schliiter regime. 
Predictions of the radial electron heat fluxes profiles according to the TFT 
and the neoclassical theory. 



is the most general coordinate transformation, under the weak-field approxi- 



mation. Moreover, Eq.(lOl) ensures that ifh^ u (X) is a solution of Eq.(94), 



so will also be h'(X) given by Eq.(lOl). h fll/ (X) and the transformed field 



h'(X) vanish as the system approaches equilibrium and, at the boundary 



const. This property is called the 



gauge invariance of the field equations. Notice that in the generalized polar 
coordinates (f , 9{, </>) with i = 1 ■ ■ ■ n — 2, r = rRo and 



X 



X 1 

X 2 
X 3 

n-l 
X n 



r cos U\ 

r sin Oi cos di 

r sin 0\ sin 62 cos 6*3 

r sin 9i ■ ■ ■ r sin Q n -2 cos < 
r sin 6% ■ ■ ■ r sin Q n -2 sin < 



(104) 



Eqs (103) take the form 



r^(0, = 
de 







where 



£ 



(X) , de v (X) 
dX^ 



dX" 



+ 



(105) 



Fig. (16) illustrates the gauge invariance property. 



Now, to perform explicitly our calculations, we have to choose a particular 
gauge i.e., to decide in which particular coordinate system we want to work. 
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Figure 13: The amplification factor F ps versus the minor radius p for a 
L-mode JET plasmas, with Z = 1. As we can see, this factor may reach 
values of order 10 2 . 



It is often convenient to work in the harmonic coordinate system [29] , defined 

J3l 



A 2 X' M = (106) 

Since 

A 2 A^ = (g XK X» x ) ;K = -Tl K g XK = -V (107) 



considering g XK qx*x«- = Eq.(107) is satisfied if = 0, which, in the 
weak-field approximation, reduces to: 

L dX^~2 L dX^ (108) 

Notice that this choice is always possible by performing the following coor- 
dinate transformations 

X » _> X'" =X fl + e^(X) (109) 

with e^{X) satisfying the differential equation 

t \k 9% _ T \ K dh KV 1 T \ K dhxK nin x 

L tVM^- L dX*~2 dX^ (11Uj 



3 In literature, the harmonic coordinate system applies to the hyperbolic P.D.E., written 
in a covariant form. Here, we indicate the harmonic coordinate system as the one satisfying 
the elliptic P.D.E., written in a covariant form. 
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Figure 14: The nonlinear Pfirsch-Schluter effect: comparison between the 
radial electron fluxes profiles according to the neoclassical theory (left pic- 
ture) and the Thermodynamic Field Theory (right picture). Although dis- 
playing similar shapes, they strongly differ in magnitude. 



with boundary conditions 



de^(X) de u (X) 



dX 1 - 



dX» 



de^(X) de u (X) 



dX l 



dX^ 



B(R ) 



x=o 



(111) 



(112) 



If h^y does not satisfy Eq.(106), we can always find an appropriate h'(X) 



that does, by performing the coordinate transformation (|109 ) with e^(X) 



and h^(0) 



satisfying Eqs (IllOh- (|112|) . Indeed, if h^ v {X) satisfies Eq.(94), with bound- 
ary conditions h fJ , u (X)\ B(Ro ^ ) = c n 
deduced that 



0, from Eqs (110)-(112), is 



dX x X K 

hlu,(X)\B(Ro) 

h'^iO) = 



const. 



and 



L 



Are 



d_KJX) 
dX x 



Kx. dh' XK (X) 
2 dX v 



(113) 



(114) 



We shall now prove that, in the limit of validity of the weak-field approxima- 
tion, i.e., up to first order in e, the boundary conditions respect the principle 
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Figure 15: The nonlinear Pfirsch-Schluter effect: comparison between the 
(radial) electron heat fluxes profiles computed by the neoclassical theory 
(left picture) and the Thermodynamic Field Theory (right picture). Also 
in this case, the profiles are qualitatively similar but, in magnitude, they 
differer by a factor, which may be of order 10 2 . 



of covariance. Let us first note that putting g^ u {X) = + h^ v {X), from 
Eq.([97} we find 

KM') = KM) + 0(e 2 ) (115) 
In particular, at the leading order, in X' = 0, we easily obtain h'(X') \x'=o= 
0. Indeed, from Eq. (|95|) we have 







(116) 



h'^(X') | X '=o= h'^Xo) + 0(e 2 ) where X* + e»(X 
Therefore, at the leading order in e, we find 
h'(X') | x , =0 = h'{-e{X )) + 0(e 2 ) = h'{X) | x=0 +0(e 2 ) = + 0(e 2 



(117) 

In the new coordinate system X /fl , /i' (X') satisfies Eq.(94) and, up to 
the first order in e, remains constant, with constant equal to c^ u , on the 
hyper-sphere of radius i?o, centered at the origin of the axes X' = 0. In- 
deed, indicating with d and d' the two boundaries in the and the X ,fl 
coordinate systems respectively, we can write 

?2 




0' 



b llv X^X v = 

^ v x'»x w 



Rl + ^I'Y + S^X^e! 1 + 0(e 2 ) (118) 
In polar coordinates, boundaries d and & are defined as 

d : r = Rq 



0' 



Ro + 0(e 1/2 ) 



(119) 
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Figure 16: The gauge invariance of the field equations: if h^ u (X) is solution 
of the field equations (satisfying the boundary conditions and 7^(0) = 0), 
so will also be tensor hL u (X) = h flu (X) + e^ u (X), with ip V (X) any tensor of 
the type plotted above (continuous lines), f indicates the normalized radial 
polar coordinate (f = r/Ro). 



From Eqs (101) and (115) we find 





= M*) 


- e^X) 


= c^u = const 




d' 


a 


a 



(120) 



On the other hand, we also have: 



9 n . 



2) ' " " > Vn-2-i 



i u n-2i ' 



+ 0(e 3 / 2 ) (121) 



2) ' " " > v n-1i 



9 n . 



+ 0{e 



3/2^ 



r'=R 



where we have taken into account that tensors h' and h^u are quantities 



of the first order in e. However, in polar coordinates, Eq.(120) reads: 



ina 

hfj, u (r, 6i,02, • • • , n -2, 



r'=R 



Efiv{ r , "l, "21' ' ' > "n-2i ' 



+ 0(e 3 / 2 )= Cflu + 0(e^)) 

r'=R 



,3/2^ 



r'=R 

(122) 
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Summarizing, If we have to solve Eq.(94) in a particular coordinate system 
where h^ v {X) is constant, say c^, on the hyper-sphere of radius Rq 
and h^ u (0) = then, after a coordinate transformation (102)-(103), also 



the new tensor h'^^X 1 ) is solution of Eq.(94) and, in the limit of validity of 



the weak-field approximation (i.e., up to the first order order in e), it also 
assumes the same constant value Cn V , on the hyper-sphere having the same 
radius Rq, in the space X' 11 . Moreover, at the origin of the axes X ,fM = 0, we 
also have h'^ u (0) = 0. In other words, up to the first order, the covariance 



principle is respected. Fig. (17) illustrates the covariance principle. 




Figure 17: The covariance principle: If h^ u {X) is solution of Eq. (94) 



satisfying the boundary conditions h^X) \b(r )= c^v = const, on a hyper- 
sphere of radius Rq, centered at the origin of the axes X = 0, and h^^O) = 
0, then, after an infinitesimal coordinate transformation, also h'^^X') is 
solution of Eq. (94) and, up to the first order in e, it assumes the same 
constant value c^ u on the hyper-sphere of radius Rq centered at the origin 
of the axes X' = 0. Moreover, h'(0) = 0. 



Notice that, with the gauge choice (109)-(112), from Eqs (96), (115) and 



(117), we derive 



L 



dX ,x X ,K 
x K dh' KU (X') 



-- + 0(e 2 ) 
\x*dti XR {X>) 



dX' x 2 dX w 
KAX')\B{R. ) = ^ + O{e Z/2 ) 
h'(0) = + O(e 2 ) 



+ 0(e 2 ) 

-- const. + 0(e 3/2 ) 



(123) 
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9 Appendix: Solutions of the Field Equations 

As seen in the last section, in a convenient choice of the coordinate sys- 
tem (the gauge choice) and after having performed an orthogonal coordi- 
nate transformation, our problem reduces to finding the solutions of Laplace 
equation. In this section we shall derive the appropriate Dirichlet boundary 
conditions. The explicit solution for the two-dimensional case, can be found 
in appendix (9.1 £ 



The Laplace equation, in two or three dimensions, reads: 

A 2 V(X) = (124) 



Solution of Eq.(124) should however satisfy the following conditions: 

a) The solution should vanish at the origin of the axes: 

M°) = ( 125 ) 

b) The solution should be constant, say with value c^ Ul on the spherical 

surface (possibly, except in a set of measure zero where it can assume 
values different from c MI/ ): 

MX)| B (jfc) = W (126) 

where B(Rq) indicates a ball of radius Ro centered at the origin of the 
axes; 

c) The solution should be invariant with respect to the permutation of the 

axes x, y and z; 

d) Inside the sphere, the solution should be of class C 2 (possibly, except in 

a set of measure zero where it should be at least of class C°). 

Now, our problem consists in re-writing these four conditions in order to 
have a well-posed Dirichlet's problem. First, we proceed by analyzing the 
simple case of Laplace's equation where the solution is of class C 2 in all 
points inside the sphere and constant on the spherical boundary, i.e.,: 

A 2 V(X)=0 where h„ v (X)\ Bm = e> + (127) 

As known, from the mean value theorem |33j . the solution is 

h^ v (X) = c^v = const. (128) 



lr The solution in a 3-dimensional thermodynamic space, is presented in ref.[9]. 
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in all points inside the sphere. Solution (128) does not satisfy the condition 
a) and for this it should be discarded. Thus the solution can be found by 
performing domain decomposition and then apply the appropriate boundary 
conditions in such a way as to obtain a harmonic function satisfying the 
conditions a), b), c) and d). However, these cuts cannot be arbitrary but 
they have: 

i) to respect the symmetry of the problem; 

ii) to be performed in such a way that, after having obtained the solution 

in a portion of the sphere, we shall then be able to reconstruct the 
entire solution, valid for the whole sphere, making use of Schwartz's 
principle (see ref 



Our solution, corresponding to the function which satisfies Eq.( 124 ) together 
with conditions a), b), c) and d), will be obtained by performing the mini- 
mum number of cuts satisfying conditions i) and ii). In this purpose, let us 
cut the sphere with a plane passing through the origin and solve the P.D.E. 



(124) with the following boundary conditions: 



u (t> a a\ f c M ,v = const, if < 9 < § . . 

V(Ao,M) = | iff <9 <l (129) 

where v^ifi, 4>) are functions of the variables 9 and <p. Note that cutting the 
sphere with a plane, passing through the origin of the sphere, and not with a 
generic surface, is the only possibility that we have if we want to respect the 
condition i). Of course without loss of generality, by performing a rotation 
of the axes, it is always possible to r each the plane z = 0. Now, den oting 
with hfj, v {r,6,<f>) the solution of Eq. (124) with boundary conditions (129) 



in the upper hemisphere < 9 < f , we are able to write the entire solution 
valid in all points of the sphere thanks to Schwartz's principle: 

I -hfj,u{r,9,(f>) iff <0<7r 



Solution (1 130|) is of class C° on the plane z = if we have 



h^(r, |,<£) = -V,(r, -,4>) i.e., h^(r, |, <f>) = (131) 



In conclusion, the solution (130), obtained by applying Schwartz's principle, 



is of class C° inside the sphere if it vanishes on the plane z = 0. This 
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constraint respects the condition a). On the other hand, to find a har- 
monic function, which vanishes on the plane z = and is constant on the 
upper spherical surface, corresponds to the solution of a well-posed Dirich- 
let's problem. In this case, the solution exists, it is unique and can be easily 



found by solving, for example, the P.D.E. (124) with the following boundary 
conditions: 

Using Schwartz's principle, the harmonic solution, having the same value 
constant on the entire spherical surface (except in a set of measure zero), 
reads 



h^(r,9,(j))= (133) 

oo 

^ J2[Pn-l(0) - P n+1 (0) + P n+ l(-l) - P n _l(-l)]P n (cOS0)( 



?1=1 

for < 9 < § , and 



r 

Rq 



K v {r,e,4>)= (134) 



oo / \ i 

V[P n -l(0) - P„+l(0) + P„+l(-l) - P n -l(-l)]Pn(cOS0) — 

^1 v^o; 



for | < # < 7r. In the last expressions, P n (cos(#)) indicate Legendre's 
polynomials. As we can see, expressions (133) and ( |134 ) do not represent 



the solution we are looking for because they satisfy the conditions a), b) 
and d) but not condition c). However, a solution, which is also invariant 
with respect to the permutation of the three axes, can be evaluated by 
cutting the sphere with three perpendicular planes, passing through the 
origin, and requiring that the solution vanishes also on the planes x = 
and y = (in addition to the plane 2 = 0). For this we have to solve 
the well-posed Dirichlet's problem in the first 3-dimensional quadrant with 
boundary conditions zero on the three planes x = 0, y = and z = 
and, constant on the spherical boundary. The entire solution, valid for all 
quadrants, will be reconstructed using Schwartz's principle. Therefore we 
seek a harmonic solution, which is constant on each spherical surface of the 
3-dimensional quadrants, with value alternatively equal to +c fJiU and to — Cuu 
(going counter clock-wise starting from the first 3-dimensional quadrant) i.e.: 



if < (f) < f 
-c p , v |ff<*<* (135) 

C Ml „ if -7T <(/><-§ 



if -5 < 6 < 
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for < 9 < § , and 



■'Mi" 



-p,,v 



if < <p < 
if f < 

if — 7T < 

if -f < 



2 

< 7T 

/»< - 

i < 



(136) 



for ^ < 9 < tt. The solution of our problem, satisfying the conditions a), b), 
c) and d), can be then obtained by applying Schwartz's principle to the ex- 
pression valid for the first 3-dimensional quadrant. The validity of the above 
method is not limited only to the three dimensional case. Indeed, it can be 
easily generalized to find the solution in a space having n-dimensions. In 



particular, for the 2-dimensional case, we have to solve Eq.(124) submitted 
to the following boundary conditions: 



if < 9 < 



if f < 
if 
if 



2 

< 7T 



-7T < < - 

< 9 < 



(137) 



and then apply Schwartz's principle. This procedure is shown in the follow- 
ing subsection. The 3-dimensional solution can be found in the appendix of 
the second part of this work |9 . 



9.1 The Two-Dimensional Problem 

In this section we shall obtain the harmonic solution, in two-dimensions, 
satisfying the conditions a), b), c) and d) mentioned in the previous section. 
The solution of this problem can be found in ref. [I]. Here, we shall re- 
obtain the same solution using the method illustrated in the last section. 



(see Fig. 18): 



Let us then start solving Eq.(124) with the following boundary conditions 



if < 9 < f 



if f < 



if 
if 



< 7T 
7T < 9 < - 
§ < 9 < 



(138) 



The solution of Eq. (124) with boundary conditions (138) can be written in 
general as [33]: 



h^vij-,9) 



Oo 

2 



+ (a„ cos n9 + b n sin n9) ( ^— ) 
7^1 V^o/ 



(139) 
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Figure 18: The 2-dimensional boundary conditions. 



with 



u nu (R o ,0) d6 

r 

Ufj,u(Ro, 6) cos[nff) d8 

K 

r 

u flv (Ro, 0) sin(n0) d9 



for n > 1 
for n > 1 



We can easily compute the integrals in Eq.(140) and obtain 

a n = (n = 0,1,2- • 

8cos(^)sin(^) 2 
K = -c uv — (n = l,2---j 



Therefore solution u a v (r, 0) can be written as 



8c^^ cos(^)sin(^) 2 
7r ^— -' n 

n=l 



sin(n#) 



which can be brought into the form 



4c uu ^ sin(2(n - 1)6) ( r y 
vr ^ 2n - 1 V #o 

n=l x 



However, solution ( 143 ) can be compacted obtaining 

"2f 2 sin(20) 



2c 



arctan 



7T 



1 -f 4 



where r 



(140) 



(141) 



(142) 



(143) 



(144) 
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or, in coordinate x and y: 



arctan 



7T 



AR 2 xy 



K - (* 2 + y 2 ) 2 



(145) 



Solution (145) is valid only in the quadrants xy > 0. The general solution, 



valid in all quadrants, can be obtained using Schwartz's principle: 



h^ v (x,y) 



2c 



TT 



arctan 



4flg 



xy 



R o 



{x 2 + y 2 ) 2 



(146) 



As seen in appendix (|8]), the coefficients c^ v are determined imposing that 
solution (146) satisfies the gauge invariance. In section ^ we obtained 



xLav where % 15 a parameter independent of the thermodynamic 
forces and L^ u indicates the Onsager matrix. In the P-S regime, parameter 
X ps is determined by imposing the supplementary condition e) reported in 
section the losses profiles are continuous functions. This condition, 
imposes that Xps = — 1- In the classical regime, this parameter can be 
determined by a simple kinetic model illustrated in ref. [llj. We obtain 
Xd = 1- 



10 Acknowledgments 

One of us (GS) would like to pay tribute to the memory of Prof. I Prigogine, 
a pioneer researcher in this area, and to Prof. R. Balescu, who gave him the 
opportunity to exchange most interesting views in different areas of plasma 
physics. 

GS would like to thank his hierarchy at the European Commission and the 
members of the EURATOM Fusion Association, Belgian State at U.L.B. 
GS is also very grateful to Dr U. Finzi, of the European Commission, for 
his continuing encouragement, for reading some sections and making helpful 
suggestions. 



References 

[1] Sonnino G. II Nuovo Cimento, 115 B, (2000), 1057. 

Sonnino G. 2001 Thermodynamic Field Theory (An Approach to 
Thermodynamics of Irreversible Processes) proceedings of the 9th 
International Workshop on Instabilities and Nonequilibrium Struc- 
tures, Vina del Mar (Chile). 



49 



Sonnino G. 2002 A Field Theory Approach to Thermodynamics of Ir- 
reversible Processes, (These d'Habilitation a Diriger des Recherches 
- H.D.R.) - Institut Non Lineaire de Nice (France). 
Sonnino G. Nuovo Cimento 118 B (2003) 1115. 
Sonnino G. Int. J. of Quantum Chemistry 98 (2004) 191. 

[2] Sonnino G. and Evslin J., Int. J. of Quantum Chemistry 107 (2007) 
968. 

Sonnino G. and Evslin J., Physics Letters A, 365 (2007) 364. 

[3] De Groot S.R. and Mazur P. 1984 Non- Equilibrium Thermodynam- 
ics, (Dover Publications, Inc., New York). 

[4] Peeters Ph. and Sonnino G. II Nuovo Cimento 115 B (2000) 1083. 

[5] Sonnino G., Nuovo Cimento, 118 B (2003) 1175. 

[6] Sonnino G. and Peeters Ph. Chaos 14 (2004) 910. 

[7] Sonnino G., Chaos 16 (2006) 023115. 

[8] Sonnino G., 33rd EPS Conference on Plasma Phys. 301 (2006) P- 
1.185. 

[9] Sonnino G. and Peeters Ph., 2007, Nonlinear Transport Processes 
in Tokamak Plasmas. Part II: The Low-collisional Regimes, to be 
submitted for publication. 

[10] Balescu R. 1988 Transport Processes in Plasmas. Vol 1. Classical 
Transport, Elsevier Science Publishers B. V. North-Holland. 

[11] Sonnino G., Physicalia, 29 no 4 (2008) 161. 

[12] Balescu R. 1988 Transport Processes in Plasmas. Vol 2. Neoclassical 
Transport, Elsevier Science Publishers B. V. North-Holland. 

[13] Prigogine I. 1947 Etude Thermodynamique des Phenomenes 
Irreversibles, (Desoer, Liege). 

[14] Prigogine I. 1954 Thermodynamics of Irreversible Processes, (John 
Wiley & Sons). 

[15] Jou D., Casas- Vazquez J. and Lebon G. 2001 Extended Irreversible 
Thermodynamics, (Springer- Verlag, Berlin Heidelberg). 



50 



Balescu R. 2005 Aspects of Anomalous Transport in Plasmas, Series 
in Plasma Physics IoP, Institute of Physics. 

Hinton F.L. and R.D. Hazeltine R.D. Rev. Mod. Phys., 48 (1976) 
239. 

Bikerton R.J., Connor and J.W. Taylor J.W. Nature, Phys. Sci. 229, 

(1971) 110. 

Kadomtsev B.B. and V.D. Shafranov V.D. Nucl. Fusion Suppl. 209 

(1972) 110. 

Wesson J. 2004 Tokamaks, (Oxford University Press). 

Hirshman S.P. and Sigmar D.J. Nuclear Fusion 24 (1981) 1079. 

Glansdorff P. and Prigogine I. 1971 Thermodynamics of Structures, 
Stability and Fluctuations, (John Wiley & Sons). 

Nishikawa K. and Wakatani M. 1999 Plasma Physics. Basic Theory 
with Fusion Applications, (Springer, Springer Series on Atoms + 
Plasmas). 

Pfirsch D. and Schluter A. 1962 Report MPI/PA/7/62, Max Planck 
Institute fhir Physik und Astrophysics, Garching. 

Haase R. 1969 Thermodynamics of Irreversible Processes, (Dover 
Publications, Inc., New York). 

Sonnino G. 2007 Plasma as Thermodynamic System, Annual report 
of the Euratom-Belgian State Fusion Association/ULB-PART. 

Sozzi C., Minardi E., Lazzaro E., Mantica P. and JET EFDA Con- 
tributors JET Profiles Analysis Based on SME, 30th EPS Conf. on 
Contr. Fusion and Plasma Phys., St. Petersburg, 27 A, (2003), 1. 

Klages R. 2003 Microscopic Chaos and Transport in Thermostated 



Dynamical Systems arXiv:nlin. CD/0309069 vl 26 Sep 2003. 



Weinberg S. 1972 Gravitation and Cosmology: Principles and Ap- 
plications of the General Theory of Relativity, (John Wiley & Sons) 

Andrews G.E., Askey R. and Roy R. 1999 Special Functions, Ency- 
clopedia of Mathematics and its Applications 71 (Cambridge Uni- 
versity Press). 



51 



[31] Courant R. and Hilbert D. 1937 Method of Mathematical Physics, 
Volume 1 Wiley Classics Editions Published in 1989, (John Wiley 
& Sons). 

[32] Morse P.M. and Feshbach H. 1953 Methods of Theoretical Physics, 
(Mc Graw-Hill Book, Inc.). 

[33] Zachmanoglou E.C. and Thoe Dale W. 1976 Introduction to Partial 
Differential Equations with Applications, (Dover, Inc., New York). 

[34] Gradshteyn I.S. and Ryzhik I.M. 1994 Table of Integrals, Series and 
Products, A. Jeffrey, Editor, (Academic Press, New York). 



52 



